

NASA Technical Memorandum 88923 
AIAA-87-0283 


Two- and Three-Dimensional Viscous 
Computations of a Hypersonic Inlet Flow 


William G. Kunik and Thomas J. Benson 
Lewis Research Center 
Cleveland, Ohio 

and 

Wing-Fai Ng and Arthur Taylor 

Virginia Polytechnic Institute and State University 

Blacksburg, Virginia 


Prepared for the 

25th Aerospace Sciences Meeting 

sponsored by the American Institute of Aeronautics and Astronautics 
Reno, Nevada, January 12-15, 1987 



=r 

St 
V ) 

I 

r- 

co 

as 


UJ 

<Tj £T> 
•H ^ 
O cm 
£3 O 

o ^ 


^r 

m 

\ 

<^» 

C 5 


-3 O 
< H!N 
as wj 
O M 

M SB 

CO 
ss 

^ u 
*3 H 

M SB 


l-J 

CJ 

CO 

u 


—4 o 
t <o 

M-l MS 
uj uj 
LB Ut 
Hi 8* 
H 33 


=5 

«S Ml. 

a 

i 

o to 

C* 25 
CH U 

m a 

tn 

— » «s 
^ pi r 
rs| a 
cn ft. 

00 *3 * 

00 o «s 

t CJ CO 
*3 «a 

H wj sb 

| CD — 

•x o 
CO CJ 3t 

CO <j 

SB H H 
' S» U, 



E- 3356 


TWO- AND THREE-DIMENSIONAL VISCOUS COMPUTATIONS 


OF A HYPERSONIC INLET FLOW 


William 6. Kunik and Thomas 3. Benson 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

and 

Wing-Fai Ng and Arthur Taylor 
Virginia Polytechnic Institute and State University 
Blacksburg, Virginia 24061 


SUMMARY 

A three-dimensional parabolized Navler-Stokes code has been used to inves- 
tigate the flow through a Mach 7.4 inlet. A two-dimensional parametric study 
of grid resolution, turbulence modeling and effect of gamma has been done and 
compared with experimental results. The results show that mesh resolution of 
the shock waves, real gas effects and turbulence length scaling are very impor- 
tant to get accurate results for hypersonic inlet flows. In addition, a three- 
dimensional calculation of the Mach 7.4 inlet has been done on a straight 
sideplate configuration. The results show that the glancing shock/boundary 
layer interaction phenomena causes significant three-dimensional flow in the 
inlet. 


INTRODUCTION 


Recent recommendations from the Aeronautical Policy Review Committee (under 
the direction of the White House Office of Science and Technology) have 
revived national Interest In hypersonic propulsion. Much of the technology 
base for this effort is expected to develop from the progressive supersonic 
and hypersonic advances in computational fluid mechanics. Several factors 
affecting the accuracy of hypersonic inlet calculations are studied in the 
current investigation by comparing the results of two-dimensional viscous cal- 
culations with experimental results. 


Hypersonic flows are characterized by high Mach numbers in the inviscid 
regions, high gradients in the vicinity of shock waves or solid surfaces, and 
thick, heated boundary layers. To accurately model the high gradients near 
shocks and along surfaces, one must employ large numbers of mesh points in the 
calculation. Most flow calculations In this regime encounter some degree of 
post-shock flow oscillations, which can be minimized either through some form 
of artificial viscosity, upwind differencing techniques, or increased grid 
resolution. To accurately model the high temperatures found in hypersonic 
flows, one needs to Include real-gas effects In the energy equation models. 
This model can take several forms from the simplest specific heats as a func- 
tion of temperature to the more complex chemical kinetic and specie models. 

To accurately model the thick boundary layers present in hypersonic flows, one 
may need to consider a variety of turbulence and transition models. In the 
present study, a supersonic parabolized Navler-Stokes solver has been used to 
investigate some of these problems. 
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Computer code modelling of the flow In supersonic and hypersonic Inlets has 
Increased recently due to Improved computer algorithms and more powerful com- 
puters such as the Cray-XMP and Cyber 205. Some researchers (refs. 1 to 3), 
have used the unsteady Reynolds-averaged Navler-Stokes (NS) equations marched 
In time to a steady state solution to describe the flow In the Inlet. At NASA 
Lewis, the Inlet has been analyzed with a zonal methodology (ref. 4); NS codes 
are used only In he transonic terminal shock region, while higher speed para- 
bolized Navler-Stokes (PNS) codes are used In the supersonic and subsonic dif- 
fusers. Through an extensive code verification effort (ref. 4), It has been 
found that the higher computational speed of the PNS codes provides the ana- 
lyst sufficient grid resolution to resolve flow phenomenon and that lack of 
sufficient grid resolution often leads to erroneous results (ref. 5). To 
attempt the same levels of grid resolution with an NS solver as Is available 
to a PNS solver would require orders of magnitude greater computing time due 
to the time-dependence of the NS solution and due to restrictions on the mag- 
nitude of the time Increment relative to the spatial grid size. 

An acknowledged dl ssadvantage of the PNS solver for high speed flows Is the 
failure of the solver when large regions of subsonic or reverse flow are 
encountered. The failure of a PNS analysis to accurately model large regions 
of separated flow In a mixed-compression Inlet may not be a serious problem, 
however, because such flows will trigger an Inlet unstart. In an Inlet 
unstart, a normal shock Is formed Inside the Inlet by the flow separation and 
this shock moves very quickly upstream. The oblique shock structure within 
the Inlet Is destroyed creating subsonic flow throughout the Inlet. The 
performance of the Inlet changes dramatically since the flow becomes subsonic 
through a very strong normal shock Instead of the multiple obliques and weak 
normal shock of a well-designed mixed-compression Inlet. Mixed compression 
Inlets are never designed with large regions of subsonic or reverse flow 
because of the Inlet unstart problem. It Is necessary In the analysis of 
mixed compression Inlets to accurately predict the occurence of reverse flows 
which trigger unstarts, but It Is not necessary to calculate the details of 
the recirculation. 

The three-dimensional supersonic viscous marching analysis used In this 
study solves the PNS equations for supersonic flow by an efficient linearized 
block Implicit scheme (ref. 6). The code has been extensively verified at 
NASA Lewis, with particular emphasis on the calculation of the glancing 
shock/boundary layer Interaction (GSBLI) (ref. 7). These studies have demon- 
strated the code's ability to properly model the complex three-dimensional phe- 
nomena present In this Interaction. It also uncovered the Importance of, and 
established the levels of grid resolution necessary to model this Interaction. 
Besides demonstrating the code's ability to match the experimental data, these 
studies have Increased our understanding of the physical processes present In 
these Interactions and their effects on Inlet performance. This Interaction 
Is very Important In the three-dimensional calculations presented here. 

TWO-DIMENSIONAL RESULTS 

Several calculations have been made of the flow through a Mach 7.4 Inlet 
and the results of the calculations have been compared to the experimental 
results. The hypersonic Inlet to be calculated was originally tested at NASA 
Ames and the experimental results are documented In reference 8. A schematic 
drawing of one configuration of the Inlet Is shown In figure 1. At Mach 7.4 
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the nominal 6.5° wedge produces a Mach 6.0 flow at the entrance to the Inlet. 

A shock Is generated by the cowl lip and this shock traverses the Internal flow 
passage and Is reflected by the ramp surface. In the vicinity of the shock 
reflection, the flowfleld was surveyed with a traversing pitot pressure probe 
to obtain pressure profiles throughout the Inlet. By varying the grid resolu- 
tion, turbulence models, and ratio of specific heats In the calculations, the 
sensitivity to these factors can be assessed. These answers can then aid 
future calculations In this speed regime and beyond. 

Figures 2 to 4 show some results from calculations of the Inlet. In each 
figure a schematic of the Inlet Is shown at the top and a comparison of the 
calculated and experimental pitot pressure profiles Is shown at the bottom. 

The location of the profiles are noted on the upper schematic. In all cases 
the results of the analysis are given by solid or dashed lines, while the 
experimental results are noted by the circles. In all of the calculations, 
the free stream conditions were set to the experimental tunnel conditions of 
Mach number equal to 7.4, total pressure equal 4.14xl0 & N/m 2 and total 
temperature equal to 811 K. All of the calculations were run fully turbu- 
lent, although there was some Indication In the test that transition occurs on 
both ramp and cowl surfaces. 

Figure 2 shows a comparison of two calculations with the experimental pitot 
profile. Considering the profile at the left, the lower part of the figure 
shows the thick boundary layer which has grown on the ramp surface; the thick- 
ness Is Indicated by the location where the curve turns vertical. Continuing 
up the curve, the sharp turn to the right Is Indicative of the cowl shock loca- 
tion while the gradual rise to the right and up Is caused by the distributed 
compression on the Inside of the cowl. Near the top of the curve, It turns 
quickly back to the left which Indicates the thickness of the boundary layer 
on the cowl. In general, the comparison between analysis and experiment Is 
excellent; the shock location and strength has been correctly modelled, as has 
the thickness of the ramp boundary layer. The boundary layer on the cowl Is 
thinner than the calculated results because of laminar to turbulent transi- 
tioning effects. 

The two calculations of figure 2 have been run with different numbers of 
mesh points In the direction normal to the ramp and cowl surfaces. When one 
hundred mesh point are used, the flow experiences severe post-shock pressure 
oscillations which are most evident In the middle profile dashed curve. Other 
researchers have experienced these oscillations In high speed flow calculations 
particularly when central differencing techniques are used. While others are 
attempting to overcome these numerical problems with special differencing 
schemes, or various types of artificial damping, this study Indicates that the 
problem can be overcome with Increased mesh resolution. This Is shown by the 
smooth solid curves on figure 2, which used three hundred mesh points In the 
normal direction. Adequate control of the post-shock oscillations was obtained 
using a minimum of two hundred radial mesh points. The resulting mesh of two 
hundred thousand mesh points took approximately three minutes of CPU time on 
the Lewis CRAY X-MP. 

Figure 3 shows the results of two calculations employing different length 
scales In the turbulence model. The dashed line shows the results using a 
McDonald-Camarata model which has been successfully used In the PNS code for 
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many different types of problems in the low supersonic speed regime (refs. 5, 

7, and 9). For the hypersonic calculations, however, this model breaks down. 
Bushnell and Beckwith found similar problems with this model in their earlier 
research (ref. 10), and using a length scale model based on this work, the 
results appear as the solid line. The comparison with experiment is much 
better and the Improved model is recommended for calculation of hypersonic 
flows. This calculation indicates possible difficulties in the calculation of 
high speed flows for which no comparison with experimental data is available; 
a turbulence model which works quite well at lower speeds gives erroneous 
results for higher Mach number flows. 

Figure 4 shows two calculations of the Mach 7.4 inlet with different values 
of the ratio of specific heats. The dashed line was calculated at gamma equal 
to 1.4, while the solid line was calculated at 1.38. A rather small change in 
this parameter has produced a large change in the pitot profiles. At the 
higher value of gamma, the shock from the cowl has moved upstream and become 
stronger. This figure indicates two important points; first, real gas effects 
which can vary the local value of gamma will have first order effects on the 
accurate calculation of the flow in high speed inlets, and second, when testing 
Inlets In facilities, it will be necessary to match flight Mach number, temper- 
ature and gamma to obtain meaningful results. Testing hypersonic Inlets in 
rocket exhaust nozzles, for example, may match flight Mach number and temper- 
ature conditions, yet fall to match gamma because of the combustion by-products 
present. This calculation shows that small differences in the value of gamma 
can lead to large differences in the assessment of inlet performance 

THREE-DIMENSIONAL RESULTS 

Once the two-dimensional parametric study was completed, it was decided 
that the inlet should be examined three-dimensional ly to determine the effect 
of the sideplates on the Internal flowfield. Previous calculations of a Mach 
5.0 inlet (ref. 11), and benchmark studies of GSBLI (refs. 7 and 9), indicate 
strong three-dimensional effects occur in rectangular inlet with sideplates. 

The experiment did not provide sufficient instrumentation for detailed study 
of three-dimensional effects. Figure 5(a) shows a schematic of the probe loca- 
tions in the spanwise direction of the inlet. In the Internal passage of the 
inlet the surveys were taken close to the centerline and therefore major three- 
dimensional effects would not have been detected. The axial locations are shown 
in figure 4. Some qualitative comparison with the throat rake data will be 
given. 

The three-dimensional calculation employed a computational mesh of one 
hundred mesh points in the Y direction and fifty points in the Z direction 
(the coordinate axes are shown In fig. 5(b)). Although the two-dimensional 
results indicated that more points are needed to control the post shock oscil- 
lations, only the overall effect of the sideplates on the two-dimensional ly of 
the flow was to be examined. To simplify the calculations, the sideplate in 
the calculation begins at the cowl lip and is perpendicular to the ramp sur- 
face. In the experiment, the sideplate was cut back as shown in figure 5(b); 
the calculation does not model this cut back. 

Figures 6 to 10 shown Mach number contours and secondary velocity vectors 
at various Y-Z planes In the Internal portion of the inlet. The bottom of the 
figures show the flowfield in a cross section of the inlet; the ramp surface 
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Is on the bottom, the cowl surface Is on the top, and the sldeplates are on both 
sides. Because of flow symmetry, only half of the Inlet was calculated. The 
side view of the Inlet Is shown at the top of each figure Indicating with a 
vertical line the relative location of the plane In the Inlet. The axial loca- 
tion given with the plot, X/HC, Is referenced to the Inlet capture height. A 
scale In shown on the left hand side for the Mach number contours and on the 
right hand side a reference vector of one-tenth of the free stream velocity Is 
shown. Boundary layers are noted by the concentration of Mach contours near a 
solid surface and shock waves are shown as a concentration of Mach contours 
away from a solid surface. Shock waves are also noted In the secondary veloc- 
ity vectors when the direction of the vectors change abruptly. 

In Figure 6 the contours are shown just upstream of the cowl lip. The 
Mach number contours show that the centerbody shock wave has passed out of the 
flow field, and that a thick boundary layer has developed on the centerbody 
surface. The secondary velocity vectors Indicate that the flow Is uniform In 
the Z-dlrectlon. Oil flow tests In reference 8 show that the flow at the cowl 
lip had been turned nearly six degrees In the Z-dlrectlon because of edge 
effects on the experimental wedge. The analysis has assumed an Infinite wedge 
and does not predict this outflow. Figure 7 shows a cross section between the 
cowl lip and throat where the cowl shock wave Is moving toward the centerbody 
surface. The figure shows that a vortex Is forming on the back side of the 
shock wave as It glances across the sldeplate boundary layer; this Inter- 
action was verified In reference 7. The vortex Is pumping the low energy fluid 
In the cowl boundary layer down along the sldeplate to the ramp surface. The 
Mach number contours show this by the thickening of the sldeplate boundary 
layer on the back side of the shock wave. Figure 7 also shows that a weak 
shock Is moving out away from the sldeplate In the Mach number contours, as 
noted by the blip In the contours on the cowl surface and In the shock wave. 
This Is developed from the boundary layer growth along the sldeplate. The 
shock wave near the sldeplate Is skewed toward the ramp surface as Indicated 
In both plots. 

Figure 8 shows the Mach number and secondary velocity plots at the point 
where the cowl shock wave reflects with the ramp surface. The vortex formed 
by the GSBLI has strengthened and Is now pumping fluid along the ramp surface 
toward the centerline of the Inlet. The boundary layer In the lower corner 
has thickened much more than In the rest of the cross section. The effect of 
the strong vortex Is seen In the boundary layer along the ramp surface by Its 
thickening to about the 50 percent point In the spanwlse direction. At this 
point the three-dimensional effects of the sldeplates have not reached the 
survey plane point which is at approximately the 30 percent point from the 
centerline. At this point in the design the cowl shock wave should have can- 
celed on the ramp surface. With the thick boundary layers on all the surfaces 
and the skewing of the shock wave near the sldeplate, the shock wave does not 
cancel at this point. The contours of figure 9 are at the point where the 
shock wave has start to reflect on the cowl surface. The vortex shown In the 
secondary velocity vectors Indicate that It Is still very strong. The vortex 
Is still pumping low energy fluid Into the lower corner. The boundary layer 
shown by the Mach number contours is very distorted In the lower corner and 
along the ramp surface. The ramp boundary layer thins toward the centerline 
of the Inlet. The secondary velocity vectors show crossflow velocities near 
the Internal passage survey plane. Figure 10 shows the contours near the Inlet 
exit. At this point the flow field Is highly distorted and the shock wave Is 
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bent across the width of the cross section. The vortex due to the GSBLI Is 
situated In the lower corner and It Is still very strong. The pumping effect 
of the vortex Is being felt over approximately 75 percent of the cross section. 

Figure 11 shows the static pressure contours along the cowl and the center- 
body surfaces. Both plots are from the sldeplate, which are on the upper part 
of both plots, to the centerline of the Inlet, which are on the lower part of 
both plots. Along the centerbody surface the cowl shock wave Impingement Is 
shown by the first concentration of contour lines. At this point the shock Is 
two-dimensional over about 60 percent of the width of the Inlet, but at the 
next shock Impingement on the centerbody the shock wave Is skewed over almost 
the entire width of the Inlet. Also, after both shock reflections on the 
centerbody an overpressure region Is formed away from the centerline as was 
seen In the GSBLI case of reference 9, The static pressure contours on the 
cowl surface show that the sldeplate Is generating a weak shock wave which Is 
distorting the distributive compression on the front part of the cowl. The 
centerbody shock Impingement shows the same effects as on the centerbody sur-. 
face. A large portion of the flow field at this point Is being affected by 
the Interaction of the cowl shock and sldeplate boundary layer. 

In figure 12 contour plots of static pressure are shown In X-Y planes, 
at selected spanwlse locations. The four planes are at; (1) the centerline 
(Z » 0); (2) the Internal survey plane (Z = 0.3); (3) the 50 percent span poln 
(Z = 0.5); and (4) the 75 percent span point (Z = 0.75). The centerline plane 
shows the same shock wave pattern as the two-dimensional results. At the 
Internal passage survey plane the shock pattern Is very similar to the center- 
line shock pattern. The only major difference Is near the last shock reflec- 
tion on the centerbody, where the effect of the sldeplate causes the shock 
Impingement point to move upstream. At the 50 percent span plane the progres- 
sion of the shock pattern upstream Is seen after the first reflection on the 
centerbody. The shock wave generated by the cowl lip does not show much change 
from the centerline shock. The final plane shows the first centerbody shock 
reflection point moves upstream of the reflection point at 50 percent. The 
reflected shock Is more highly Inclined to the ramp and strikes the cowl 
further upstream. This change In the shock wave pattern with spanwlse location 
may help explain a phenomenon noted In the experimental report. A Schllerln 
photograph of this region shows the cowl shock striking the centerbody and two 
waves reflecting back towards the cowl. The computational results Indicate 
that this may be due to the curvature of the shock wave seen In figure 12. 

This does not affect the cowl shock wave because the three-dimensional ly of 
the flow field has not developed enough at this point and It Intensifies after 
the shock wave boundary layer Interaction on the ramp surface. 

The three-dimensional flow features seen In the previous figures cause the 
performance of this Inlet to be less than predicted two-dimensional ly. Fig- 
ure 13 shows pitot pressure contours at the exit of the Inlet as calculated by 
the three-dimensional analysis. Only half of the Inlet Is shown In this fig- 
ure, with the sidewall at the left. The GSBLI has caused the pitot contours 
to be highly distorted In the corner formed by the ramp and sidewall. The 
effects of the Interaction are seen to extend over more than 50 percent of the 
cross section. The Instrumentation of the experiment was concentrated along 
the centerline except near the throat. Figure 14 shows a comparison of the 
analysis and experimental results In this region. The comparison can only be 
made qualitatively because of Insufficient Information In the experimental 
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report. However, the figure shows that both analysis and experimental Indi- 
cate strong variations In pitot pressure with span; and that both Indicate 
higher pitot pressures near the cowl surface near the sidewall. 

CONCLUSIONS 

From the two-dimensional results presented here three major conclusions 
can made. First, to adequately model the shock wave structure and to elimi- 
nate post shock oscillations In high speed flow fields a large number of mesh 
points are required. The results of this study Indicate that two hundred 
points In the Y-dlrectlon are adequate to control post shock oscillations In 
this analysis. The Increased number of mesh points can still be calculate very 
quickly with the analysis. Other researchers have Indicated that the same 
results can be obtained by using alternate differencing techniques. Second, 
to model t.he boundary layer growth In these flows, which are very large, proper 
scaling of the turbulence Is Important. Very good agreement with the experi- 
mental data Is obtained using a mixing length model with proper corrections 
for high speed flow. Third, real gas effects will become very Important In 
accurately predicting the correct pressure distributions In high speed, high 
temperature flows. The results presented here show that changing gamma over 
the entire flow field dramatically change the free stream pressure levels. A 
more sophisticated model for real gas effects should be Incorporated later. 

The three-dimensional calculation presented In this report have led to the 
following conclusions. The glancing shock/boundary layer Interaction produces 
a strong vortex that persists throughout the Internal passage of the Inlet. 

This vortex redistributes the boundary layer Into the lower corner of the Inlet 
and then proceeds to distort the ramp boundary layer. To achieve a canceled 
shock wave In the Inlet, the effect of the boundary layers on the shock wave 
Interaction at solid boundaries must be taken into account. Also, the effect 
of the sldeplates on the shock wave necessitates that Inlet design must be 
based on an understanding of the three-dimensional nature of the flow 
field. 


REFERENCES 

1. White, M.E.; Drummond, J.P.; and Kumar, A.: Evolution and Status of CFD 

Techniques for Scramjet Applications. AIAA Paper 86-0160, Jan. 1986. 

2. Knight, D.D.: Improved Numerical Simulation of High Speed Inlets Using 

the Navler-Stokes Equations. AIAA Paper 80-0383, Jan. 1980. 

3. Campbell, A.F.; and Forester, C.K.: Evaluation of a Method for Analyzing 

the Aperture Region of Two-Dimensional External Compression Inlets. AIAA 
Paper 85-3072, Oct. 1985. 

4. Anderson, B.H.: Three-Dimensional Viscous Design Methodology for Advanced 

Technology Aircraft Supersonic Inlet Systems. AIAA Paper 84-0194, Jan. 
1984. 

5. Benson, T.J.; and Anderson, B.H.: Validation of a Three-Dimensional 

Viscous Analysis of Axlsymmetrlc Supersonic Inlet Flow Fields. AIAA Paper 
83-0137, Jan. 1983. 


7 


6. Buggeln, R.C., et al.: Development of a Three-Dimensional Supersonic 

Inlet Flow Analysis. NASA CR-3218, 1980. 

7. Anderson, B.H.; and Benson, T.J.: Numerical Solution to the Glancing 

Sidewall Oblique Shock Wave/Turbulent Boundary Layer Interaction In Three 
Dimension. AIAA Paper 83-0136, Jan. 11583. 

' ^ 

8. Gnos,A.V., et al.: Investigation of Flow Fields within Large-Scale 

Hypersonic Inlet Models. NASA TN D-7150, 1973. 

9. Benson, T.J.; and Anderson, B.H.: A Study of Three Dimensional Shock Wave 

Turbulent Boundary Layer Interactions Present Within Aircraft Inlets. 

AIAA Paper 84-1558, June 1984. 

10. Bushnel 1 ,D.M. ; and Beckwith, I.E.: Calculation of Nonequlllbrlum 

Hypersonic Turbulent Boundary Layers and Comparison with Experimental 
Data. AIAA J., Vol . 8, No. 8, Aug. 1970, pp. 1462-1469. 

11. Benson, T.J.: Three-Dimensional Viscous Calculation of Flow In a Mach 

5.0 Hypersonic Inlet. AIAA Paper 86-1461, June 1986. 


8 



NORMAL DISTANCE FROM 
HALL Y. CM 


ORIGINAL PAGE IS 
Of POOR QUALITY 



FIGURE 1.- MACH 7.R HYPERSONIC INLET GEOMETRY. 


NUMBER OF NORMAL MESH POINTS 
300 
100 



^ 

r \ 



A 

A 

j' 

// 

I 

/' 

. . . . 



PITOT PRESSURE RATIO. PS/PSREF 


C0-M-Z20U 


FIGURE 2. - EFFECT OF GRID RESOLUTION ON INLET CALCULATION. 
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